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Abstract 

Inspired by soliton models, we propose a description of static particles in terms of 
Riemannian 4-manifolds with self-dual Weyl tensor. For electrically charged particles, 
the 4-manifolds are non-compact and asymptotically fibred by circles over physical 3- 
space. This is akin to the Kaluza-Klein description of electromagnetism, except that we 
exchange the roles of magnetic and electric fields, and only assume the bundle structure 
asymptotically, away from the core of the particle in question. We identify the Chern 
class of the circle bundle at infinity with minus the electric charge and the signature of 
the 4-manifold with the baryon number. Electrically neutral particles are described by 
compact 4-manifolds. We illustrate our approach by studying the Taub-NUT manifold 
as a model for the electron, the Atiyah-Hitchin manifold as a model for the proton, 
CP 2 with the Fubini-Study metric as a model for the neutron, and S 4 with its standard 
metric as a model for the neutrino. 

1 Introduction 

Geometry and the quantum mechanics of particles have an uneasy relationship which is 
why general relativity is hard to incorporate into quantum field theory. String theory is an 



ambitious and remarkable attempt at unification with many successes but the final theory 
has proved mysterious and elusive. 

Einstein and Bohr fought a long battle on this front and Bohr was generally deemed to have 
won, with the Copenhagen interpretation of quantum mechanics accepted. But Einstein's 
belief in the role of geometry made a partial come-back with the adoption of gauge theories 
as models of particle physics. 

A more modest and limited role for geometry in nuclear physics was proposed by Skyrme 
[1] with the solitonic model of baryons, i.e. proton, neutron and nuclei, now known as 
Skyrmions. These have been shown to be approximate models of the physical baryons 
occurring in gauge theories of quarks and gluons, and have been extensively studied [2] with 
considerable success. 

In this paper we explore a geometric model of particles which is inspired by Skyrme's idea 
but with potential applications to both baryonic and leptonic particle physics. Our model 
differs from the Skyrme model in that it uses Riemannian geometry rather than field theory, 
and so it is closer in spirit to Einstein's ideas. Another key difference is that we absorb the 
Kaluza-Klein idea of an extra circle dimension to incorporate electromagnetism. However, 
we exchange the roles of electricity and magnetism relative to the standard Kaluza-Klein 
approach, so that the extra circle dimension is magnetic rather than electric, and it is the 
electric charge that is topologically quantised by the famous Dirac argument. We will also, 
initially, ignore time and dynamics, focusing on purely static models. 

Our geometric models will therefore be 4-dimensional Riemannian manifolds. We require 
the manifolds to be oriented and complete but not usually compact: we use non-compact 
manifolds to model electrically charged particles, and compact manifolds for neutral particles. 
Both baryon number and electric charge will be encoded in the topology, with baryon number 
(at least provisionally) identified with the signature of the 4-manifold 1 and electric charge 
with minus the Chern class of an asymptotic fibration by circles. In particular, the number 
of protons and the number of neutrons will therefore be determined topologically. 

The manifolds will also have an 'asymptotic' structure which captures their relation to 
physical 3-space. The non-compact manifolds we consider have an asymptotic region which 
is fibred over physical 3-space, so no additional structure is required. For the compact 
(neutral) models, however, we fix a distinguished embedded surface X where the 4-manifold 
M intersects physical 3-space. We also call M\X the inside of M. 

For single particles, the symmetry group of rotations should fix all of the above data. It 
should act isometrically on the 4-manifold and preserve the asymptotic structure. In the non- 
compact cases this means that it should be a bundle map in the asymptotic region, covering 
the usual S'0(3)-action on physical 3-space. In the compact cases it should preserve the 
distinguished surface X. In order to capture the fermionic nature of the particles considered 

1 The definition of signature for non-compact manifolds is reviewed in Sect. 8.1. 
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in this paper we also require spin structures on the non-compact manifolds and on the inside 
(in the sense defined above) of compact manifolds. Moreover, the lift of the rotation group 
action to the spin bundle (over the entire manifold in the non-compact case and over the 
inside in the compact case) should necessarily be an SU (2)-action, and this is what we mean 
by saying that our models are fermionic. 

A key restriction on (the conformal classes of) our manifolds is that they are self-dual. 
Recall that the Riemann curvature is made up of the Ricci tensor plus the Weyl tensor W, 
which is conformally invariant. In dimension 4, W is the sum of self-dual and anti-self-dual 
parts 

w = w + ®w- . (1.1) 

A 4-manifold M is said to be self-dual if W~ = 0. These manifolds are precisely those 
which have twistor spaces in the sense of Penrose [3]. The latter are 3-dimensional complex 
manifolds Z with a real S^-fibration over M. The complex structure of Z (together with a 
real involution which is the antipodal map on each S 2 ) encodes the entire conformal structure 
of M. Even the Einstein metric can be captured by complex data on Z. 

A simply-connected self-dual 4-manifold, for which the Ricci tensor is also zero, is a hy- 
perkahler manifold, whose structure group reduces to SU (2). It is a complex Kahler manifold 
for an S^-family of complex structures, but for any of these complex structures, and for the 
complex orientation, it would be anti-self-dual. Since we want self-dual manifolds we choose 
the opposite orientation. 

While some of our particles, including the proton, will be modelled by hyperkahler mani- 
folds, we do not want to be so restrictive. Instead we will only require our 4-manifold models 
of particles to be self-dual and Einstein, so there can be a non-zero scalar curvature. Our 
model for the neutron is of this type, distinguishing it from the proton. The neutron will of 
course also have electric charge zero. 

We should point out that reversing the orientation of a 4-manifold turns a self-dual manifold 
into an anti-self-dual one. This should be interpreted as giving the geometric model of an 
anti-particle. The existence of anti-particles follows from CPT invariance, and our models 
are compatible with this. 

Self-dual 4-manifolds are, in many ways, the 4-dimensional analogue of Riemann surfaces, 
with H 2 replacing Hi in homology. In particular there are theorems [4] which assert that 
such manifolds admit connected sums although, unlike in the case of Riemann surfaces, there 
are restrictions on when this is possible. Such connected sums model composite objects, like 
nuclei. Although we focus at present on static particles we do envisage a deformation theory, 
using the moduli space of self-dual manifolds, which could underlie particle interactions. 

Fortunately a lot is now known about self-dual 4-manifolds with many metrics explicitly 
calculated. This makes it possible to put forward some definite models for the proton and 
neutron. Even though our ideas are inspired by Skyrme's theory of baryons, it turns out that 
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geometric models of leptons, i.e. the electron and (electron-)neutrino, are even simpler, and 
we shall describe them too in the class of self-dual manifolds. Thus, somewhat surprisingly, 
our framework of self-dual manifolds allows us to describe baryons and leptons in a unified 
fashion. 

The language and spirit of our model for particles is close to that of general relativity 
and suggests the possibility of a unification with gravity, but we do not address this issue 
here. In particular, we do not specify an action functional. Instead, we focus on how our 
model describes general features of particles such as their various quantum numbers. We are 
aware that a description of elementary particles as 4-dimensional Riemannian manifolds is 
radically different from established treatments in terms of quantum field theory. What we 
aim to show in this paper is that such a geometric approach is possible, and that it has some 
surprising and attractive features, such as the possibility of describing the electron and the 
proton in one framework. While we do propose definite identifications of certain 4-manifolds 
with specific particles in Sects. 3-5 of this paper, these should be seen as illustrations of the 
geometric approach, not necessarily as final proposals. 

The paper is organised as follows. In Sect. 2 we outline the genesis of our geometric 
models of particles, starting with the Skyrme model of baryons. Electrically charged particles 
are necessarily described by non-compact 4-manifolds in our approach, and in Sect. 3 we 
explain how to model the electron and proton in terms of the Taub-NUT and Atiyah-Hitchin 
manifolds, respectively. Neutral particles are described by compact 4-manifolds, and this is 
discussed in Sects. 4-5. We propose CP 2 as a model for the neutron and S 4 as a model for 
the neutrino. These are the simplest choices, but we also discuss some more sophisticated 
versions. In Sect. 6 we describe how our particle models glue into empty space, and how 
the particles may interact with each other. Sect. 7 contains an outline of how our geometric 
models capture the spinorial nature of the particles they describe. In Sect. 8 we give the 
dictionary which translates topological properties of 4-manifolds into the electric charge and 
baryon number of particles, and discuss in some detail how these charges are related to fields 
and densities used in conventional Lagrangian models of particle physics. Sect. 9 contains 
our conclusion and some ideas for follow-up work. Conventions and calculations are collected 
in appendices A,B and C. 



2 From Skyrmions to 4-manifolds 

We begin by spelling out in detail how the Skyrme model suggests our 4-manifold model. 
The Skyrme model is based on a group- valued field from M 3 , 

U:R 3 ->G, (2.1) 

where the Lie group G is usually taken to be SU(2), and U(x) — > 1 as \x\ — > oo. The degree 
of U as a map S 3 — > SU (2) is identified with baryon number. The minima of the Skyrme 
energy, for each baryon number, are called Skyrmions. 



4 



Skyrmions are free to rotate, both in physical space and through conjugation by elements 
of SU(2). Quantising this motion gives the Skyrmions spin and electric charge. The proton 
and neutron, for example, are distinct quantum states of the essentially unique Skyrmion of 
degree 1. 

In [5] it was shown how to generate such Skyrme fields naturally by starting with an SU (2) 
Yang-Mills gauge field on IR 4 and calculating the holonomy along the 4th direction. Suitable 
asymptotic behaviour on IR 4 guarantees a well-defined map U. Although this construction 
does not preserve the respective energy functionals it does provide a good way of using 
instantons on IR 4 (i.e. self-dual gauge fields) to construct approximate minima of the Skyrme 
energy. It also identifies instanton number with the Skyrme degree. See also the recent 
papers [6, 7] where the difference between the Yang-Mills and Skyrme energy functionals is 
interpreted as due to an infinite tower of mesons. 

Since the Yang-Mills energy functional in dimension 4 is conformally invariant we could 
replace the decomposition 

]R 4 = K 3 x IR 1 (2.2) 

by 

S A \ S 2 = H 3 x S 1 , (2.3) 

where H 3 is hyperbolic 3-space. In fact we can vary the curvature of H 3 provided we rescale 
the circle S l the opposite way, so that large circles correspond to almost flat H 3 . We can 
now fix a gauge field on S A \ S 2 and take the holonomy round the circles. There are some 
technicalities (due to base-points) which we shall ignore but basically we expect to end up 
with a Skyrmion on H 3 , an idea which has been explored in [8]. 

Now replace H 3 x S* 1 by any Riemannian 4-manifold M which is asymptotically fibred 
by circles over IR 3 . This is the kind of Kaluza-Klein 4-manifold we are going to consider. 
An SU(2) gauge field on M would then give a Skyrme field on 'the quotient of M by S 1 \ 
Since we do not want to assume there is a global circle fibration, this Skyrme field will only 
be defined asymptotically outside some 'core'. But an oriented 4-manifold has two natural 
£77(2) bundles over it, the two spin bundles £ + and S~ (assuming M is a spin-manifold, i.e. 
W2(M) = 0). Picking one of these, say £ + , we then get from the connection on S + a natural 
construction of an asymptotic Skyrme field on IR 3 . 

This is roughly the genesis of our idea to model particles by 4-manifolds, but the topology of 
these asymptotic Skyrme fields does not quite fit, and would not give integer baryon numbers 
as defined in our model. The reason lies in a fundamental difference between the topology 
of gauge bundles which can have arbitrary instanton number (or second Chern class) and 
the topology of tangent bundles, where there are divisibility theorems. For example the 
first Pontrjagin class of a 4-manifold (compact and oriented) is divisible by 3 and then gives 
the signature. An example is CP 2 which has signature 1 and Pontrjagin class 3 (times the 
generator of H 4 {CF 2 )). 
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In the Skyrme model the basic idea is that baryon number is identified with the degree 
of the map U in (2.1), or equivalently with the instanton number (or second Chern class) 
of the SU{2) bundle over IR 4 . This differs from the 4-manifold model we want to explore, 
where baryon number is identified with the signature of the 4-manifold. The signature is 
additive under taking connected sums of 4-manifolds [9], and this captures the additivity of 
baryon number for composites of particles, for example, fusion of nuclei. The integrality of 
the signature is linked to it being an index of an elliptic operator. This means we are in the 
realm of K-theory rather than cohomology. A consequence of this change of viewpoint is 
that the geometry of the 4-manifold model is important for us, but we will not try to define 
a global 3-dimensional Skyrme field U. 

Recall that in the Skyrme model, baryon number is cohomological and electric charge arises 
at the quantum level. For our 4-manifold model, electric charge is cohomological, arising, 
as already explained, from the first Chern class of the asymptotic S^-fibration, while baryon 
number as just indicated should be seen as an index. 

To sum up our discussion, we see that our model goes beyond the Skyrme model in aiming 
to understand topologically both the basic integer physical invariants, baryon number and 
electric charge. The two models are different, but possibly dual in a suitable sense. We hope 
to explore this in detail at a later stage. 

3 Models for the electron and proton 

Models for the basic particles should exhibit a high degree of symmetry and we expect the 
rotation group SO (3) of IR 3 , or its double cover £77(2), to act as isometries. For electrically 
charged particles, we take our geometric models to be non-compact hyperkahler manifolds. 
We also assume that the volume grows with the third power of the radius, to allow for an 
interpretation of the asymptotic region in terms of physical 3-space. As recently shown in 
[10], this forces the hyperkahler manifold to be ALF. We are therefore looking for rotationally 
symmetric and complete ALF hyperkahler manifolds. There are just two possibilities: 

1. The Taub-NUT manifold [11, 12] depending on a positive parameter m (interpreted as 
mass in the gravitational context). For brevity we denote it by TN. 

2. The Atiyah-Hitchin manifold, the (simply-connected double cover of the) moduli space 
of centred £?7(2)-monopoles of charge two [13]. For brevity we denote it by AH. 

Note that we could also single out TN and AH among non-compact, complete and rota- 
tionally symmetric hyperkahler manifolds by demanding that the £77(2)- (or £0(3)-) action 
rotates the complex structures, see our discussion following (B.12). This turns out to play a 
role in recovering the usual rotation action on physical 3-space in the asymptotic region of 
our geometric models, as discussed in Sect. 7. 
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Both TN and AH can be parametrised in terms of a radial coordinate r and angular 
coordinates on SU(2) (for TN) or SO(3)/7*2 (for AH). Details are given in Appendix B. In 
terms of the left-invariant 1-forms defined in (B.4), the metrics of both TN and AH can be 
written as 

ds 2 = f(r) 2 dr 2 + a(r) 2 r? 2 + b(r) 2 r] 2 2 + c(r) 2 r] 2 3 , (3.1) 
with the functions /, a, b, c satisfying the self-duality equations 

y^ = (&-c) 2 -o 2 , + cyclic, (3.2) 

where + cyclic means we add the two further equations obtained by cyclic permutation of 
a, b, c. We adopt the convention 

/W = -^, (3-3) 

where (for reasons that will emerge later) the radial coordinate r has the range [0, oo) for 
TN and [n, oo) for AH. The self-duality equations (3.2) become 

^ = -l(a 2 -(6-c) 2 ), 
dr 2rc y y 11 

£_-L ((? _(._»«. (3.4, 
This system has solutions in terms of elementary functions 



a(r) — b(r) — r \ e -\ and c— — . m , (3.5) 

y j k > y r v^+F 

with parameters e,m > 0, associated to the TN manifold. The topology is that of M 4 , and 
as e — > the metric tends to the flat metric. For e > the manifold is asymptotic to an S 1 
fib re-bundle over IR 3 with the length of the circle being Airm/ \fe. There is a U(l) symmetry 
acting along the fibres, with just one fixed point at the origin, r = 0. The whole isometry 
group is U(2). As e — > the [/(l)-action becomes the scalar action on C 2 . The complex 
orientation of C 2 determines the orientation of TN as a self-dual manifold; this is opposite 
to the orientation given by any of the complex structures in the hyperkahler family, see 
Appendix A. 6 for a discussion. At infinity, the [/(l)-action gives the standard Hopf line 
bundle over CP 1 with Chern class +1; details are given in Appendix A. 

The TN metric with coefficient functions (3.5) has the following behaviour under scaling 
by non- vanishing real numbers a, j3: 

r— >— r, m-^a[3m, e — > a 2 e , then ds 2 — > [3 2 ds 2 . (3.6) 
a 



7 



We use rescaling by a to set e = 1, and rescaling by f3 to set m = 2 from now on. This 
amounts to picking a unit of length for the radial coordinate r and to fixing an overall scale 
for the metric. Our choice is motivated by the asymptotic form of the AH metric, to be 
discussed below. Note that, with this choice, the length of the asymptotic circle, in the 
length units chosen, is 8n. 

The solution which gives rise to AH has the asymptotic form, for large r, 



2 2 
a{r) ~ b(r) ~ r\l 1 and c(r) ~ . (3.7) 



r /i_ 



r 



These asymptotic expressions, a TN metric with e = l,m = —2, also satisfy (3.4). However 
a(r) is not actually equal to b(r), and r only extends down to n. For r near n, (3.7) is a 
poor approximation. Instead, one finds the leading terms 

1 1 

a{r) ~ 2(r — 7r) , 6(r) ~ ir H — (r — n) , c(r) ~ — n H — (r — n) , (3.8) 

which we will need later in this paper. 

The manifold AH is the complement of MP 2 (the real projective plane) embedded in CP 2 , 
and the complex orientation of CP 2 determines the orientation of AH as a self-dual manifold. 
As for TN this is opposite to the orientation given by any of the complex structures in the 
hyperkahler family; see our discussion in Appendix A. 6. AH has an 5*0(3) symmetry with 
just one 2-dimensional orbit at r = ir, which is a minimal 2-sphere. We refer to this minimal 
2-sphere, which is the totally imaginary conic in CP 2 and determined by z\ + z\ + z\ = 
in the homogeneous coordinates introduced in Sect. 4.1, as the core. Asymptotically, the 
manifold is fibred by circles. As further discussed below, neither the circles nor the base 
space of this asymptotic fibration are oriented because of a ^-identification, given explicitly 
in (B.7). 

The manifold TN is usually interpreted as the geometry of a Dirac monopole at the origin 
[14, 15, 16]. For us, with electric and magnetic charges reversed, it has to be interpreted as 
an electrically charged particle. Since the signature of TN is zero (we discuss this further 
in Sect. 8) the particle is leptonic. We therefore interpret TN as a model for the electron. 
Down on R 3 , after factoring by U(l), any 2-sphere surrounding the origin has an electric 
flux emerging from it due to the electron, which carries charge —1. This implies that there 
is a sign change in going from the Chern class to the electric charge. 

The manifold AH has the opposite asymptotic behaviour with a sign change for m and 
an orientation change (see Appendix A) and so would lead us naturally to expect electric 
charge +1. Also, the topology at the core is different, with a 2-sphere instead of a point. As 
a result, AH has signature 1 (again discussed further in Sect. 8) and looks like the model we 
want for the proton (rather than the positron). 
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However things are not quite that simple, as we shall now explain. The 'asymptotic 
boundary' of AH is not S 3 as for TN but the boundary of a tubular neighbourhood of IRP 2 
in CP 2 , which is S* 3 divided by a cyclic group of order 4. Moreover the base of this unoriented 
circle fibration is MP 2 , not S* 2 , and is non-orientable. This means that the 3-manifold which 
is the base of the asymptotic fibration is not M 3 , and has a fundamental group of order 2. It 
is not orientable. 

This might seem to be a disaster, but we shall argue that, while unexpected, it is not as 
bad as it looks. The most convincing argument in its favour is to show that the electric 
charge is well-defined and equal to +1 as hoped. This is done in detail in Appendix A. 
The lack of orientability in physical 3-space should be thought of as follows. The lack of 
orientation in IR 3 locally is compensated by a corresponding ambiguity in the sign of electric 
charge (non-orientability of the circle fibres). Physically the geometric orientation is not felt. 

4 The neutron 

4.1 Complex projective plane 

Having put forward a definite proposal for the proton we now have to face the neutron. 
Since the neutron has no electric charge any non-compact model would need to have a 
trivial asymptotic circle fibration. The 4-manifold should have signature 1 and it should 
resemble the AH model of the proton in its SO (3) orbit structure. However, the latter 
requirement rules out asymptotically trivial circle bundles over physical 3-space since the 
generic 5*0(3) orbits would be 2-dimensional in that case. We therefore consider compact 
4-manifolds. In fact, there is an obvious choice which is just the complex projective plane 
CP 2 with its Fubini-Study metric (and its natural complex orientation). This is a self-dual 
manifold of positive scalar curvature. 

CP 2 has even more symmetry than we need, since it is acted on by SU(3). The rotation 
group SO (3) sits inside as the subgroup that preserves the real structure given by complex 
conjugation. This preserves MP 2 as in the case of the proton, so we fix this IRP 2 as the 
distinguished surface where the 4-manifold intersects physical 3-space, thus breaking the 
symmetry to S0(3) . The global SU(3) symmetry might give us a link to quarks, but this 
remains to be explored. 

The Fubini-Study metric is often written in coordinates which exhibit the invariance under 
U(2) C SU(3). This brings out the parallels with TN [17] but is not the symmetry we want 
in our neutron model. For the interpretation of CP 2 as a neutron and for a comparison 
with the AH model of the proton, we need to write the Fubini-Study metric in coordinates 
adapted to the S'0(3)-action, which is discussed in [18] and [19]. We write the results of [18] 
in the conventions used in our discussion of the AH metric. 

In terms of homogeneous coordinates z e C 3 \{0} (with the identification z ~ Xz, A G C*) 
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the Fubini-Study metric on CP 2 is 



(41) 

For calculations we can fix \z\ = 1 and parametrise 

z = e ia Rz , (4.2) 

where R G 5*0(3) can, in turn, be parametrised in terms of Euler angles as shown in Appendix 
B.l. The reference vector Zq (which depends on one parameter) should be a unit vector and 
we can assume, by adjusting the phase a if necessary, that its real and imaginary parts are 
orthogonal. For our purpose, it is convenient to single out the 1-axis and pick 




z = I a 2 I , a\ + a\ = 1 . (4.3) 
\ia 3 J 

Then we find the following expression for the Fubini-Study metric in terms of the left- 
invariant forms (B.4) on SO (3) (for details of an analogous calculation we refer the reader 
to [18]): 

ds 2 = da\ + da\ + (a\ - a 2 .) 2 rj{ + a^rfe + alvl ■ ( 4 -4) 

Parametrising 

a 2 = cos(^- + -J, a 3 = sm^- + -J, pE [0, -\ , (4.5) 
we obtain, finally, the Fubini-Study metric in the form 

ds 2 = \dp 2 + sin 2 prji + sin 2 (| + |) V 2 2 + cos 2 (J + ^) % 2 . (4.6) 



There is a simple interpretation of the geometry of CP 2 and its orbit structure in terms of 
orientated ellipses up to scale [20], which is useful for comparison with the AH metric. As 
already exploited above, we can adjust the phase in the homogeneous coordinate z = u + iv 
(no longer fixed to satisfy \z\ = 1) so that the real vectors u and v are orthogonal: if 
z 2 = this is automatic, and if z 2 ^ we multiply by a unit complex number to set 
Im(z 2 ) = 2u ■ v = and Re(z 2 ) = u 2 — v 2 < (we pick the negative sign to agree with the 
choice made in (4.5) above). We can interpret v and u as the major- and minor-axis of an 
ellipse. This ellipse is only determined up to scale (we can still rescale z by any positive real 
number) but it is orientated. The totally degenerate case z = is excluded by the definition 
of homogeneous coordinates, but circles [z 2 — or |u| = |i>|) and lines = 0) can occur. 

In terms of our parametrisation (4.2), the reference vector (4.3) and the definition of p in 
(4.5), we see that for p = the ellipse is a circle and for p — | it degenerates to a line. For 
generic values of p, the S'0(3)-orbit is SO(3)/7j 2 , with the Z 2 generated by the 180°-rotation 
about the 1-axis, but for p = the orbit is a 2-sphere and for p — | the orbit is MP 2 . This 
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is the same as the orbit structure of AH compactified by an IRP 2 at infinity, although the 
metric is of course different. 

We also note that the Kahler form 



\z 2 > 


dz l A dz 




z^dz A z t dz 




z\ 


4 



(4.7) 



takes the simple form 

uj = cos p r/2 A rjs — sin p dp A rji , (4.8) 

which should be compared to the expression given in [17] in coordinates adapted to the 
U(2) symmetry of CP 2 . The form to is invariant under the 180°-rotation about the 1-axis 
and hence well-defined on the generic S'0(3)-orbits. It is manifestly closed, but not exact: 
for p > we can write u = d(cospr]i) but this expression is not valid on the exceptional 
S'0(3)-orbit where p = 0, since rji is not well-defined there. 

The Kahler form is self-dual with respect to the complex orientation (with volume element 
dV = u 2 = — sin(2p) dp A 7]i A i]2 A i] 3 ). Since it is closed, it is also harmonic. The existence 
of a non-exact harmonic, self-dual form on CP 2 follows from the fact that the signature of 
CP 2 is 1. In view of our interpretation of signature as baryonic charge one might expect 
there to be a baryonic interpretation of oj. We return to this question when discussing the 
AH model of the proton further in Sect. 8. 

4.2 Hitchin's one-parameter family of Einstein metrics 

If the CP 2 -model for the neutron turns out to be too naive, there is a more sophisticated 
variant which could be explored. This arises from the sequence of self-dual Einstein manifolds 
M(k), for k > 2, studied by Hitchin [21]. The manifolds M(k) for even k > 4 are all defined 
on the same space as our proton model, namely CP 2 but with WP 2 removed. M(4) is CP 2 
with the Fubini-Study metric and all the metrics on M(k), k > 4 and even, are incomplete 
on the open set in CP 2 \ IRP 2 , but can be completed to metrics on CP 2 with a conical 
singularity of angle An/(k — 2) along MP 2 . For odd k > 3 the manifolds are defined on S* 4 , 
with M(3) being S* 4 with its standard metric. This time the metrics of M(k), k > 3 and 
odd, are incomplete on the open set in ,S 4 \MP 2 , but can be completed to metrics on S* 4 with 
a conical singularity of angle 2n/{k — 2) along WP 2 . 

The sequence of conical manifolds M(k) for even k > 4 and starting with CP 2 , has de- 
creasing scalar curvature and converges to AH as k — > 00. It may turn out that some other 
value of k gives a better model for the neutron than k = 4. Note that for k > 4 the conical 
singularity breaks the symmetry down to 5*0(3). Even for k = 4 we shall see later that other 
factors break the symmetry in this way. 

Hitchin also pointed out [21, 22] that the family M(k) can be extended to real parameter 
values. For any real k > 0, M{k) is related to the moduli space of centred SU (2) monopoles 
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over hyperbolic space of curvature — l/p 2 , where p = (k — 4)/4. When k is not an integer, 
the conical angle is not a rational multiple of ix and M(k) is not an orbifold. Consequently, 
the explicit methods of [21] do not then apply. Nonetheless, having A; as a real parameter 
gives useful room for manoeuvre in modelling the neutron and may provide contact with 
conventional nuclear models. In particular, 1/k may play a role as a small parameter that 
controls the breaking of isospin symmetry. 

The forthcoming paper [23] contains a signature formula for Riemannian manifolds with 
conical singularities, like the Hitchin manifolds M(k). We summarise that result in Sect. 8.2. 

5 The neutrino 

Having put forward models of the electron, proton and neutron, it is then natural to look for 
a similar model of the neutrino. Since it has no electric charge it should, like the neutron, 
be modelled by a compact manifold. It should have symmetry similar to the U (2) symmetry 
of the electron and should have positive curvature. It should have zero baryon number, that 
is, vanishing signature. 

Just as CP 2 is the most obvious model for the neutron, the standard 4-sphere, 5 4 , is the 
most obvious model for the neutrino. Again this has more symmetry than we need, 5*0(5) 
instead of 5*0(3). Just as a distinguished MP 2 in CP 2 picks out the smaller 50(3) symmetry, 
so a distinguished 5 2 in 5 4 is needed to cut down the symmetry of the neutrino. 

To exhibit the symmetry, we parametrise 5 4 in terms of vectors x G M? and y G R 2 
satisfying the constraint 

x-x + y-y = 1 . (5.1) 

The metric is then 

ds 2 = dx-dx + dy-dy . (5.2) 

The group 50(3) x 50(2) acts in the obvious way on the pair of vectors (x,y) G 5 4 and 
preserves the metric. In order to compare with other metrics discussed in this paper we 
parametrise 

x = sin p (sin 6 cos <j), sin 6 sin 0, cos 6) , 

y = cos p (cos x, sin x) , (5.3) 

with p G [0, |], x G [0, 2n) and the usual range for the polar coordinates 9, <f> on 5 2 , and find 
the expression 

ds 2 = dp 2 + sin 2 p (rjl + r/ 2 ) + cos 2 p d\ 2 ■ (5.4) 

Here we used that d6 2 + sin 2 9 d(j) 2 = i] 2 + i] 2 in terms of the left-invariant 1-forms defined in 
(B.4). The generic 50(3) x 50(2) orbit is 5 2 x 5 1 , but this collapses to 5 1 when p = and 
to 5 2 when p = n/2. 
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Note that S* 4 is conformally flat, so the Weyl tensor vanishes and is trivially self-dual. S* 4 
has no middle-dimensional homology, so the signature is zero, and hence our model neutrino 
has zero baryon number, as required. Since S A also has an orientation-reversing isometry, 
our model seems to suggest that the neutrino coincides with the anti-neutrino. For this and 
other reasons (see later) our choice of S* 4 is very tentative and provisional. As with the CP 2 
model for the neutron it should be regarded at present as a prototype. 

To address the symmetry breaking issue, and several others, we will in the next section 
discuss how our various models are supposed to fit into conventional 3-space. 

6 Attaching the models to space 

So far our models are abstract objects, 4-manifolds on their own, which are supposed to 
model four basic particles of nature. How are we to view them in the real world? 

Let us begin with the easiest case, that of the electron. Thought of originally as the 
Dirac monopole, the idea is well-known. We consider Kaluza-Klein space as a Riemannian 
4-manifold with a circle action. Away from matter this space is assumed to be a circle bundle 
over R 3 . Outside a given region in IR 3 which is electrically neutral the bundle is assumed to 
be (topologically) just the product space. If a region is electrically charged the circle bundle 
over the boundary is supposed to have a Chern class equal to minus the charge. If we start 
from the vacuum, then inserting one electron amounts to attaching a truncated version of 
TN to the boundary. This truncation turns the idealized model into a more realistic model of 
a particle. If other particles are present, TN will be an approximation to the precise metric. 
This approximation is some measure of the force exerted on the electron. In a dynamic 
theory, forces should emerge from the equations, a task for the future. 

Next let us move on to the proton. This is similar, using a truncated version of AH. 
However, as pointed out earlier, the circle fibration is now not oriented so that the asymptotic 
3-space is not IR 3 (unless our region contains equal numbers of protons and anti-protons). 

The model for the neutron is compact, so there is no way to attach it to a boundary. Instead 
we propose that our model neutron (a copy of CP 2 ) intersects our 4-space in a surface. This 
surface should project to a surface in 3-space which is the 'boundary' of the neutron as seen 
by an observer. Since we want to keep the neutron similar to the proton (except for the 
charge) it seems reasonable to take this surface to be a copy ofRP 2 . But since the charge 
is now zero the circle bundle over the surface should be trivial. This could arise as follows. 
Pick a point in M 3 and blow it up to give an MP 2 , so that we are modifying 3-space in this 
neighbourhood. Keep the circle bundle trivial, so making the charge zero. Lift this MP 2 into 
the total space of the circle bundle and let the CP 2 neutron model intersect 4-space in this 
IRP 2 , which is the distinguished MP 2 in CP 2 (hence breaking the larger symmetry of CP 2 ). In 
this construction the metric on CP 2 need not be changed. The only change that is needed is 
the change in metric on the background 4-space got from the blowing up process in 3-space. 
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Finally we come to the neutrino. From the other examples it is clear what is required. 
This time we just take a 3-ball in IR 3 with boundary a 2-sphere. The circle bundle over it is 
trivial and we lift the sphere to the total 4-space. We now require the S A neutrino model to 
intersect our 4-space in the chosen 2-sphere. The choice of the 2-sphere in S* 4 again breaks 
the symmetry, down to 5*0(3) x SO(2). The 2-sphere need not be a great (geodesic) sphere 
and this provides a parameter to play with. 

Note that in this case we could reinterpret the picture as the surgery that kills off the circle 
and leaves a 2-sphere. This means that for the electron, the proton and the neutrino we can 
still think of our 'space' as a 4-manifold. But this does not seem to work for the neutron 
where we have to settle for a 4-space with intersecting components like a complex algebraic 
surface with double curves. 

The four models that we have proposed for the four basic particles should be geometrically 
related in some way to account for the process of beta decay in which a neutron breaks up 
into a proton, an electron and an anti-neutrino. The opposite asymptotic behaviour of AH 
and TN is a good start but the difference in the asymptotic fundamental groups presents a 
problem. This suggests that the model of the neutrino should somehow bridge the gap and 
it argues against the simplicity of the 4-sphere. We hope to pursue this question. 

7 Spin 1/2 

In all our models we have a natural action of the symmetry group of rotations (SO (3) or 
SU(2)) preserving the metric and the 'asymptotics', the details of which differ according to 
the cases. For the neutral, compact models we interpret 'asymptotic' to mean the behaviour 
near the distinguished surface where the 4-manifold intersects 3-space, which is either an 
IRP 2 or an S 2 . For the electrically charged, non-compact models we have an asymptotic 
fibration by circles over physical 3-space; rotations preserve this fibration and induce an 
5'0(3)-action on the base. The hyperkahler structures on TN [24] and AH [25] can be 
used to construct Cartesian coordinates on physical 3-space with the physically required 
transformation properties under spatial rotations. Here we make essential use of the fact 
that the complex structures on TN and AH transform as vectors under rotations, as explained 
at the beginning of Sect. 3. 

For a model to represent a particle of spin 1/2 we must include the data necessary to lift the 
rotation group action to an SU (2) action, and to construct its 2-dimensional representation. 
To achieve this for non-compact (electrically charged) models we require a spin structure 
on the 4-manifold while for compact (neutral) models we only require a spin structure on 
the inside, obtained by removing a distinguished surface X from the 4-manifold. In this 
short section we explain why, in each of the models considered in this paper, the lift of the 
rotation group action to the spin bundle is an SU (2)-action. We do not attempt to construct 
naturally associated spin 1/2 representations, but comment on how this may be done. 
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For the TN model of the electron, there is nothing to do since the rotation group action 
on TN is an SU (2)-action. For the neutron model we view the required data as the com- 
pactification of the proton model AH, not just topologically, but also with the action of the 
rotation group. In particular the MP 2 at infinity is part of the data. We now simply require 
the extra data of a spin structure on the inside CP 2 \ MP 2 , i.e., on AH. 

It might be thought that, since the spin structure on AH is unique, there is nothing gained 
by the additional data, but this is to ignore the interaction with the symmetries. We must 
now require the rotation group to lift to the spin bundle, and this may require us to pass to 
SU(2), in which case we label the model as fermionic. Otherwise we call it bosonic. To see 
that, in principle, either case could occur, consider the two manifolds Y 1 = 50(3) x R and 
Y 2 = M 3 x R. We take the left-translation action of SO (3) on Y\ and the standard action on 
Y 2 . For Fx, the tangent bundle is trivial and we can choose the trivial spin bundle (though 
since Y 1 is not simply-connected there is another choice). The 50(3)-action extends without 
going to SU(2), so Y x is bosonic in our terminology. For Y 2 the fixed point at the origin of 
IR 3 means that we can only lift to the spin bundle after passing to SU (2), so Y 2 is fermionic. 

For AH we have to show that, with its action of 50(3), it is fermionic. There are several 
ways to do this. Perhaps the simplest (in line with the example Y\ above) is to note that the 
action is not free and that the isotropy group of a point on the core is 5*0(2). To lift even 
this subgroup to the spin bundle over the fixed point requires us to go to the double cover. 

A comparison between our models for the electron and proton is illuminating. As we 
pointed out, SU(2) acts on the electron but only SO (3) acts on the proton. Thus the 
fermionic natures of our two models differ. In one case it is inherent in the symmetry while 
in the other it is geometric or topological. 

Our tentative model for the neutrino is just the round 4-sphere, with a distinguished 2- 
sphere at infinity given by the decomposition IR 5 = IR 3 x M 2 and the corresponding action of 
5*0(3). The inside, got by removing S 2 , has an infinite cyclic fundamental group so there 
are infinitely many spin structures, but only one extends to the whole of S 4 . If we pick this, 
then it is easy to see that the lift of SO (3) to the spin structure requires us to pass to SU (2) 
(for example we can use an 50(2) fixing a point at infinity and argue as with the proton). 
Thus our model of the neutrino is fermionic. 

Our discussion so far shows that our geometric approach furnishes fermionic models, but 
it does not establish that they necessarily give spin 1/2. This requires constructing the 2- 
dimensional representation of 577(2) and relating it to the asymptotic region. We expect that 
the required 2-dimensional representations can be constructed in terms of eigenspaces of the 
Dirac operator on our model manifolds, possibly twisted by a [/(l)-bundle with curvature 
proportional to one of the harmonic 2-forms discussed in Sects. 4.1 and 8. 
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8 Charges, energies and fluxes 



8.1 General remarks 

So far we have focussed on topological and geometrical features of our models and explained 
how they describe general properties of particles - like baryon number, electric charge and 
location in space. We want to keep an open mind about how our geometric models make 
quantitative contact with the physics of elementary particles. In particular, we do not assume 
that this should necessarily happen in the standard framework of Lagrangian field theory, 
where dynamics, conservation laws and even the quantum theory are all derived from an 
action functional. 

The purpose of this section is to illustrate that our geometric models for particles nev- 
ertheless contain natural candidates for the kind of quantities which arise in Lagrangian 
models, like energy density and electric fields. We show that electric charge as defined in 
our model can be represented by a harmonic 2-form, thus making contact with the usual 
description of electric flux. One important feature of the densities and fields considered in 
this section is that they are defined on the 4-manifold so that they can only be interpreted 
as conventional spatial densities and fields in the asymptotic region of the 4-manifold which 
canonically projects down to physical 3-space. 

We begin with a summary of the topological quantities and their physical interpretation for 
each of the 4-manifolds considered thus far. For compact manifolds the electric charge is zero 
and for non-compact manifolds it is minus the self-intersection number X 2 of the manifold 
X representing infinity in their compactification. The compactification is CP 2 for both TN 
and AH, but X = CP 1 for TN with self-intersection number X 2 = 1, while X = MP 2 for 
AH, with self-intersection number X 2 = — 1; see Appendix A. 6 for details. 

The baryon number is identified with the signature of the 4-manifold. For a non-compact 
oriented manifold, the signature is defined as the signature of the image of the compactly 
supported cohomology in the full cohomology [26] . The topology of TN is that of IR 4 so the 
signature vanishes. The signature of AH is 1. This follows from the fact that its 2-dimensional 
homology is generated by the core 2-sphere, and that the self-intersection number of this 
2-sphere is positive (in fact equal to +4). The same argument applies to the sequence of 
Hitchin manifolds M(k) for k > 4 and even, reviewed in Sect. 4.2, which are all compact 
and topologically equivalent to CP 2 . 

In Table 1 we list the electric charge and baryon number as well as the Euler characteristic 
X and the squared L 2 -norm ||-R|| 2 of the Riemann curvature for the four 4-manifolds mainly 
discussed in this paper. The Euler characteristic is homotopy invariant, so can be computed 
for TN and AH by noting that the former retracts to a point and the latter to a 2-sphere. 
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Table 1: Geometric properties of 4-manifolds and their physical interpretation 



8.2 Signature, Euler characteristic and L 2 -norm of the Riemann curvature 

For compact Riemannian 4-manifolds there are formulae for the Euler characteristic and 
signature in terms of integrals over the 4-manifold involving the Riemann curvature, which 
we shall review below. In the non-compact cases, these bulk contributions need to be sup- 
plemented by boundary integrals and (for the signature) a subtle spectral contribution (rj- 
invariant), see e.g. [12] and [23] for a summary. For manifolds with conical singularities like 
the sequence of Hitchin manifolds reviewed in Sect. 4.2, a signature formula was recently 
found [23], which we also review below. We now discuss the bulk contributions, relegating 
most detailed calculations to Appendix B.2. 

Writing R for the Riemann tensor as in Appendix B, we define the squared L 2 -norm of R 
(for compact and non-compact manifolds) as 

n^n 2 = / y / R ab A*R ab , (8.1) 

jM a<b 

where *R a b is the Hodge dual of R a b- The form that integrates to the first Pontrjagin class 
on a compact manifold is 

pi = ^J2 RabA Rab ' ( 8 - 2 ) 

a<b 

and the form which integrates to the signature r in the compact case is 

S=\pu (8-3) 



so that 



T (M)= f S=\f Pi = T^[ 52 R* A R*. (8.4) 



a<b 

The form that integrates to the Euler characteristic x m the compact case is 



^Ys^Rab^Rcd- (8.5) 



16tt 2 

a<b 
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Note that, with our conventions and for compact Einstein manifolds [27], 

x(M) = L e = ^ m? ' (8 - 6) 

which determines \ \R\\ 2 for CP 2 and S* 4 in terms of their topology. 

The Riemann curvature on a 4-manifold may also be thought of as a mapping of 2-forms. 
Exploiting the fact that the space A 2 of 2-forms on an oriented 4-manifold decomposes into 
±-eigenspaces A* 1 of the Hodge star operator *, we get a corresponding decomposition of the 
Riemann curvature into irreducible pieces [12] 

Here W ± are the self-dual and anti-self-dual parts of the Weyl tensor, s is the scalar curvature 
and B amounts to the tracefree part of the Ricci curvature. Then the signature of a compact 
manifold can also be expressed as [27] 

r(M) = I ^(||H/ + || 2 -||H/-|| 2 ). (8.8) 

For a self-dual manifold W~ = so, up to the factor 127T 2 , r is given by the L 2 -norm of the 
Weyl tensor W + , and is non-negative. 

Since the metrics on both TN and AH are hyperkahler, B and s vanish, so the full Riemann 
curvature is self-dual. Therefore the bulk contribution to both the signature and the Euler 
characteristic can be expressed in terms of the L 2 -norm of the Riemann curvature, which 
equals the L 2 -norm of W + . In Appendix B.2, we show that 

\\R\\l N = 8?r 2 , H^HIh = 16tt 2 . (8.9) 

Using the self-duality of the Riemann curvature, we deduce that the bulk contributions to 
the Euler characteristic are both in agreement with the topological results listed above: 

r 8tt 2 r i6tt 2 

L e= ^ =u L e= ^ =2 - i8M) 

The bulk contributions to the signatures, on the other hand, turn out to be fractional: 

T / = l^ = 3' J An S= ^ = 3- (8 - U) 

As shown in [23], the fall-off of the spin connection and curvature imply that the boundary 
integrals do not contribute in the limit to either the Euler characteristic or the signature. 
However, the fractional values of the bulk integrals for the signature show that there must 
be a non-zero contribution from the 77-invariant. 
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In [23], a signature formula is derived for Riemannian 4-manifolds with conical singular 
metrics. As reviewed in Sect. 4.2, the Hitchin manifolds M(k) are of this type, and M(oo) can 
be identified with AH. The general formula derived in [23] for the signature of a Riemannian 
4-manifold M whose metric has a conical singularity of angle 2tt/k along a surface X is 

tw = si* -H 1 -;?)**■ (812) 

where pi is the form (8.2) which integrates to the Pontrjagin class, and X 2 is the self- 
intersection number of X in M. 

Since the signature r is interpreted as the baryon number in our model, it is tempting to 
interpret the integrand of the bulk contribution to the signature as a baryon number density. 
Our calculations above show that this cannot be the whole story, since the signature r also 
receives a contribution from the 77-invariant of the boundary. We nevertheless compute and 
plot the integrands of the bulk contributions to the signature below. Since overall factors 
are not important here, we look at the squared L 2 -norm (8.1) of the Riemann curvature. 
In the notation of Appendix B, the integrand of (8.1) in both the TN and AH case can be 
written as 

r dF 

dF A r]i A r/2 A i] 3 = —— — dV , (8.13) 
ao z c dr 

where dV is the volume element, so that the combination 

VI < 8 - i4 > 

ao z c dr 

may be interpreted as a density. We plot this density for TN and AH in Fig. 1. It is finite 
at, respectively, the origin and the core, and for large r falls off like 1/r 6 . In the AH case we 
compute the functions a, b, c and hence F numerically; in the TN case we have the explicit 
formula 

r iF 24 (8.15) 



ab 2 c dr (r + 2f 



8.3 Electric and baryonic fluxes 

We now construct fields on TN and AH which carry the electric flux measured by the 
self-intersection numbers X 2 tabulated in Table I. We do this by systematically studying 
rotationally symmetric harmonic forms on both TN and AH, with most of the details given 
in Appendix C. 

One finds [28, 29, 30, 31] that the only square-integrable harmonic and rotationally sym- 
metric 2-form on TN is, up to an overall arbitrary constant, 

^ = Gts * A * + vvw dr ^)= d (vh ^ ■ (8 - 16) 
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Figure 1: The density (8.14) for TN (left) and AH (right) 



We claim that this form is a harmonic representative of the Poincare dual of the surface 
X = CP 1 at infinity in TN, the surface parametrised by 9 and 0. Our calculation also gives 
an alternative computation of the electric charge as minus the self-intersection number of 
CP 1 in CP 2 . Let p be this self-intersection number, then 



is Poincare dual to CP 1 since 



uj cv i= p . (8.18) 

CP 1 



However, by the definition of the self-intersection numbers in terms of cohomology, we also 
have 

%'Au CP i=p. (8.19) 



/ 



'TN 

Evaluating the integral in Appendix C we find 

p 2 =p, (8.20) 

so that p — 1, confirming that, with the self-dual orientation, the self-intersection of CP 1 in 
CP 2 is 1. 

Since the 2-form — cu CP i is harmonic and since its total flux through infinity equals the 
electric charge we interpret it as the electric field of the electron. Although we have adopted a 
viewpoint dual to standard electromagnetism, with a purely spatial 2-form being interpreted 
as electric rather than magnetic flux, it is interesting that the self-dual 2-form (8.16) also 
contains a term which allows for a conventional electric interpretation: when we contract 
u C¥ i with the vector field along the fibres (using r) 3 (d^) = 1), we obtain a purely radial 
field which, asymptotically, falls off like 1/r 2 . 
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The integral (8.19) is the squared L 2 -norm of the electric field. Ignoring overall factors 
and working with we write 

IM1| 2 = J PtdV, (8.21) 
with dV defined as in (B.17), and interpret the integrand 



(r + 2Y 



.22) 



as an electric energy density. For comparison with the AH case below, we plot the profile 
function g£(r) — r/(r + 2) and the energy density p% in Fig. 2. 




0.12 
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Figure 2: The 'electric' profile function g£ and the energy density p% 

Turning now to AH, we review in Appendix C why there are only two S'0(3)-invariant 
harmonic forms £lf on AH which respect the identification (B.7). They have the structure 



Of = Gfr l2 Ar ]3 + dGf A Vl , 



123) 



for functions Gf of r which satisfy the ordinary differential equations for gf in (C.3). Using 
the asymptotic formulae (C.13) and (C.14) one shows [24, 32, 33] that only Gf is finite at 
the core and decays at infinity. In fact, the solution decays exponentially fast at infinity, 
with the leading term proportional to e -r / 2 . We normalise (n) = 1, so that near the core 



Gi(r) ~ 1 -(r — 7r) 2 , for (r — n) small. 



(8.24) 



The 2-form Clf is not dual to the surface X 

in terms of electric flux. However, it is dual to the core CP 1 in AH. With 

V 



a 



An 



at infinity in AH and has no interpretation 

(8.25) 
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one finds 

n coic =p. (8.26) 



/ 

J CO 



As a check, we use this form to compute the self-intersection number of the core. Recalling 
that the volume of SO(3)/Z 2 is — 4ir 2 , and using (C.15), we have 



I ri corc Afi corc = -|^(-47r 2 )(0 2 -l 2 ) = ^, (8.27) 



so that 

j=P, (8-28) 

confirming p = 4 for the self-dual orientation (this agrees with [13], where the anti-self-dual 
orientation gives the opposite sign). 

The harmonic form (4. 8) on CP 2 is related to the signature of CP 2 through the fact that 
the signature of a compact 4-manifold equals the difference of the dimensions of the spaces 
of self-dual and anti-self-dual harmonic representatives of the second de Rham cohomology. 
It seems likely that the (up to scale) unique bounded and rotationally symmetric self-dual 
harmonic form f2 cor c on AH is related to the self-dual harmonic form (4.8) via the sequence 
of Hitchin manifolds reviewed in Sect. 4.2. 

Since the existence of f2 corc on AH is linked to the signature of AH being 1, and since 
signature represents baryon number in our approach, it may be possible to interpret the 
detailed structure of fi core in baryonic terms. The exponential fall-off exhibited by Q COTC is 
reminiscent of the proton's pion field in the Yukawa description of the nuclear force. In 
Fig. 3 we plot the numerically computed profile function and the associated energy 
density (C.17), which turns out to be 

Pt = 2<gi! . (8.29) 

This is finite at the core and decays exponentially according to e~ r for large r. 

To find a harmonic 2-form on AH which can play the role of the proton's electric field 
we need to go to a branched cover of AH, denoted AH. The metric on the branched cover 
is not smooth at the core, but this will not affect the following calculations near infinity. 
Dropping the requirement that forms are invariant under the identification (B.7), we find 
that the closure condition on the 2-forms (C.ll) has only one further solution which is finite 
at the core and which remains bounded for large r. This is the 2-form 

^3 = G3 A r] 2 + dG 3 A r] 3 , (8.30) 

with a radial function G 3 satisfying (C.3). The solution vanishes at the core as 

G 3 (r) ~ CVr^ , for (r - vr) small , (8.31) 



22 



1.0 - 
0.9 - 
0.8- 
0.7 - 
0.6- 
0.5- 
0.4- 
0.3 - 
0.2- 
0.1 - 



0.005 - 



10 12 14 16 



Figure 3: The 'baryonic' profile function G± and the energy density p{ 



for some constant C. The large r behaviour is 

G£(r)~C 



r-2 



(8.32) 



where C is another constant. This leads to an anti-self-dual form on AH, which is square- 
integrable (as we shall show below) and which has not been considered previously. It is our 
candidate for the electric field of the proton. 

As a manifold, AH compactifies to S 2 x S 2 , as discussed in Appendices A. 4 and A. 5. The 
surface X = WF 2 which compactifies AH to CP 2 becomes the anti-diagonal S 2 in S 2 x S 2 . 
Choosing C = 1 in (8.32), a harmonic representative of the Poincare dual class to this S 2 is 



n 



g2 



Then 



47T d 



il S 2= P . 



S 2 



and since the volume of 50(3) is — 8n 2 by (B.8), we also compute 

„2 



rt$2 a 



AH 



16tt 2 



-8tt 2 )(1 2 -0 2 



' 2 



Hence 



which is solved by p = 
intersection number of 



pr 
' 2" 



P: 



(8.33) 
(8.34) 

(8.35) 
(8.36) 



-2. Dividing by 2, to take us back to MP 2 , we confirm the self- 



m 



as —1, again for the self-dual orientation. 



Since the 2- form —^52 on AH is harmonic and since its total flux through infinity, suitably 
interpreted, equals the electric charge, we think of —^52 as representing the electric field 
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of the proton. The comments made after (8.16) about the possibility of recovering a con- 
ventional electric field by contracting the electric 2-form on TN with the vector field along 
asymptotic fibres apply suitably modified, to AH. 

In Fig. 4 we plot the numerically computed profile function G 3 for C — 1 (as defined in 
(8.31)) and the associated energy density (C.17), which turns out to be 

(8-37) 

This function has a (r — 7r) _1 singularity at r = ir, and falls off like r -4 for large r as in the 
TN case. 




Figure 4: The 'electric' profile function G 3 and the energy density p 3 



9 Conclusion and outlook 

In introducing and illustrating geometric models of matter in this paper we have concentrated 
on general, global properties of particles like baryon number and electric charge. The most 
striking phenomenological success of our geometric approach is its prediction of precisely two 
stable electrically charged particles, one leptonic and one baryonic, with opposite electric 
charges. 

An important theme throughout this paper is the implementation of rotational symmetry 
in our models. The rotation action preserves the metric and, in the compact cases, the 
distinguished surface where the 4-manifold intersects physical 3-space. We have shown that 
our models are fermionic in the sense that the lift of the rotation action to spin bundles (over 
the 'inside' of the 4-manifold in the compact cases) is necessarily an S77(2)-action, but we 
have left the explicit construction of spin 1/2 states for future work. 
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Non-vanishing electric charge and baryon number give rise to harmonic 2-forms on our 
model manifolds. These have allowed us to make contact with the conventional description 
of charged particles in terms of associated fields and fluxes. 

There are various geometrically natural candidates for measuring energy or mass in our 
models, but we have not committed ourselves to any particular energy functional in this pa- 
per. In the absence of an energy measure, we are not able to make contact with experimental 
data about particle masses or forces between particles. 

We end our summary with some observations about relative scales predicted by our models 
for charged particles. The AH model relates three scale parameters: the size of the core 
(radius ~ it in our geometric units), the size of the asymptotic circles (radius 2 in geometric 
units) and the scale 1 implicit in the exponential corrections to the asymptotic form of the 
coefficient functions (3.7) for the AH metric, which are proportional to e _r (see [34, 35] for a 
discussion of these corrections in the context of magnetic monopoles). Interpreting the core 
radius as the proton radius and the scale in the exponential decay as the Compton wavelength 
of the pion, we find that our model correctly assigns the same order of magnitude to those 
two quantities. The details are not quite right (experimentally, the proton radius is just over 
half of A^) but this is not surprising given our very rudimentary understanding of how our 
model relates to experiment. In any case, these considerations suggest that we should pick 
a length unit of about 1 fermi or 1CT 15 m in the AH model. 

Since the asymptotic fibration by circles arises for all electrically charged particles, we 
expect the size of the asymptotic circle to have a purely electric interpretation, and we also 
assume that it is the same for both AH and TN (this was implicit in the way we fixed 
scale and units for TN). One natural guess, at least for TN, is to relate the size of the 
asymptotic circle to the classical electron radius. Our models would then equate the orders 
of magnitude of the classical electron radius with the Compton wavelength of the pion, which 
is phenomenologically right. It is an attractive feature of TN as a model for the electron 
that it has a length scale, but no core structure. 

Many avenues remain to be explored. In our geometric approach, fusion and fission of 
nuclei as well as decay processes involving both baryons and leptons (like the beta decay 
of the neutron) should have a description in terms of gluing and deformation of self-dual 
4-manifolds. In order to study masses and binding energies of particles we need to pick an 
energy measure. This could involve the norms of curvatures and harmonic forms discussed 
here, but may take a less conventional form. Our ideas about spin 1/2 need to be fleshed 
out. The Dirac operator on our model manifolds is likely to play a role in combining energy 
and spin considerations, and Seiberg-Witten theory on the model manifolds seems relevant, 
too. In order to move beyond static consideration, time needs to be introduced. 

The list of open issues may seem daunting, but in each case the geometric framework 
introduced here suggests natural lines of attack. We hope to pursue them in future work. 
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Appendices 

A Geometry and sign conventions 
A.l Self-intersection numbers 

Let L be a complex line bundle over a compact oriented surface X. Then 

Cl {L)[X]=X 2 (A.l) 

where X 2 denotes the self-intersection number of the zero section in the total space of L. 
This is essentially one of the definitions of the first Chern class. Note that both X and L are 
oriented, with the fibres having the natural orientation of the complex numbers. However the 
self-intersection X 2 depends only on the orientation of the total space, since (—X) 2 = X 2 . 

A.2 The Hopf bundle 

Let us review the definition of the Hopf bundle H, the standard line bundle over the complex 
projective line, paying particular attention to orientations. We start with C 2 , and CP 1 as 
the space of complex lines through the origin. This complex line bundle is defined to be 
the dual of the line bundle H . The reason for this apparently perverse choice is that H has 
holomorphic sections (given by C 2 ) and its dual H* does not. This is clear from the exact 
sequence of holomorphic vector bundles over CP 1 : 

0^H*^C 2 ^H^0. (A.2) 

Since holomorphic intersections are always positive, the first Chern class of H gives +1 when 
evaluated on the fundamental class of CP 1 . 

If we compactify C 2 by adding a CP 1 at infinity then the self-intersection number of this 
line is clearly +1, so that the normal bundle is H . 

A. 3 Non-orientable surfaces 

Now assume that X is a non-orientable surface in an oriented 4-manifold M. This defines a 
mod 2 homology class in M and so it would appear that its self-intersection number is only 
defined modulo 2. However a more careful look shows that we can define an integer self- 
intersection. There are several (equivalent) ways to see this. First we note that the tangent 
bundle and the normal bundle of X in M have isomorphic orientation real line bundles L. 



26 



Then both the Euler class of the normal bundle and the fundamental class of X are defined 
as classes 'twisted' by L. The evaluation (A.l) therefore makes sense and defines X 2 . An 
alternative way is to pass to the double covering of a neighbourhood of X in M, where X 
acquires an orientation, compute the self-intersection there and divide by 2. A third way 
is to deform X into a transverse surface X and, near each intersection point of X and X, 
choose an orientation of X and the corresponding orientation of X given by the deformation. 
Now compute the local intersection numbers of X and X and sum them up. 

A. 4 An example 

The example we want to study is that of the real projective plane X = WP 2 inside the 
complex projective plane M = CP 2 . We give M its standard complex orientation. We plan 
to show that X 2 — — 1. 

To do this we will use the double branched covering of CP 2 given by the product of two 
copies of the complex projective line CP 1 . We can identify CP 2 with the symmetric product 
of the two copies, branched along the diagonal D. This map is holomorphic and compatible 
with the action of £X(2,C). Now introduce a metric on CP 1 identifying it with a 2-sphere, 
and consider in CP 1 x CP 1 the graph D' of the antipodal map. The image of D' in CP 2 is 
the standard embedding of MP 2 into CP 2 . Outside the branch locus D the map from the 
product of the two 2-spheres to CP 2 is a double covering so we can compute X 2 as half of 
{D') 2 . This is most easily done by observing that, in the cohomology of the product, with x 
and y the two generators of H 2 , the classes of D and D' are 

D = x + y, D' = x-y. (A.3) 

Since x 2 = y 2 = and xy — 1 it follows that 

D 2 = 2 and (D'f = -2 (A.4) 

(while D.D' = agreeing with the fact that D and D' are disjoint). Dividing by 2 to get 
back down to CP 2 we end up with X 2 = — 1 as stated. 

As a check we can also use a deformation of D' in M by choosing a different metric on 
CP 1 giving a different antipodal map. In particular, consider the 2-sphere as the boundary 
of the unit ball in IR 3 and move the origin so that antipodal points are now the end points of 
lines through the new origin. The line joining the two origins defines a unique intersection 
point of the two copies of MP 2 in CP 2 . Since this line acquires different orientations from 
the two origins the intersection number is —1. 

A third sign check is to observe that the normal bundle to the diagonal in the product 
of the two 2-spheres is the tangent bundle with Euler number 2, while for the antidiagonal 
(graph of the antipodal map) it is the cotangent bundle with Euler number —2. 
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A. 5 More on the example 

It is enlightening to analyze the SO(3) action on the spaces in the example above. In addition 
we can introduce the branched double covering CP 2 — > S 4 given by complex conjugation 
[36, 37]. This is branched over MP 2 . 

Thus we have maps 

S 2 x S 2 ->• CP 2 ->• S 4 , (A.5) 

which are compatible with the 5*0(3) actions. Here S A is best thought of as the space of real 
3x3 symmetric matrices of zero trace and fixed norm. 

The generic orbits are, in each case, 3-dimensional but there are two exceptional orbits of 
dimension 2. Inside the product S* 2 x S* 2 they are the diagonal and anti-diagonal (both S 2 ), 
in S* 4 they are both MP 2 , while in CP 2 one is S* 2 and one is MP 2 . Each of the two maps is 
branched along a surface: S 2 for the first map and IRP 2 for the second. 

To examine the self-intersection numbers of the exceptional orbits it is useful to note two 
general rules: 

A. If X in M is double covered by X' in M' then (X') 2 = 2X 2 . 

B. If X is a branch locus in the double covering then (X') 2 = \X 2 . 

Rule A is obvious and rule B is familiar in algebraic geometry, but it applies more generally 
even in the non-orientable case. 

Using rules A and B we can now check what happens in our two maps. In S A the two 
]RP 2 s have self-intersection numbers +2 and —2. In S 2 x S 2 the two 2-spheres also have self- 
intersection numbers +2 and —2. In CP 2 the 2-sphere, which is a conic, has self-intersection 
number +4, while MP 2 has self-intersection number —1. 

A. 6 Self-dual manifolds 

We want to analyze various sign conventions in relation to self-dual manifolds, and begin by 
checking our sign conventions for two examples of complex Kahler surfaces. Referring to our 
discussion of the signature and Euler characteristic in Sect. 8.2, we note the examples: 

1. The K3 surface with the Yau metric, which makes it hyperkahler. With its complex 
orientation its signature is —16 and its Euler characteristic is 24. With the opposite 
orientation, the metric is self-dual and the signature is positive. 

2. CP 2 with the Fubini-Study metric is self-dual for its complex orientation [12], agreeing 
with the positive signature +1. 
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We will also be interested in non-compact examples similar to the above. The argument 
that a hyperkahler manifold is self-dual for the orientation opposite to the complex orienta- 
tion (given by one of the family of complex Kahler structures) is purely local. It just depends 
on the fact that the bundle of self-dual 2-forms for these complex orientations is flat. More 
generally, on any Kahler manifold, the bundle of self-dual 2-forms is the direct sum of the 
canonical line bundle K, its dual and a trivial bundle generated by the Kahler form. 

The first non-compact example that interests us is the Taub-NUT manifold TN. This 
has the isometry group U(2), its topology (and symmetry) is that of C 2 with the central 
£7(1) giving the scalar action. The other U(l) actions, inside SU(2), define the complex 
hyperkahler structures which have the opposite orientation. Hence TN is self-dual for the 
orientation that becomes that of C 2 in the limit when the parameter e goes to 0, as discussed 
after (3.5). 

The second example is that of the (simply-connected version of the) Atiyah-Hitchin man- 
ifold AH, which is an open subspace of CP 2 got by removing MP 2 . As shown in [13], it has 
a complete hyperkahler metric. This is self-dual for the orientation opposite to the complex 
orientations given by the hyperkahler metric. But this is just the orientation given by its 
complex structure as an open subset of CP 2 . This can be checked directly, but it is best 
seen by using the results of Hitchin [21], reviewed in Sect. 4.2, which give a whole sequence 
of self-dual manifolds on the same space, starting from the Fubini-Study metric of CP 2 and 
converging to AH. 

In each case, for TN and AH, we have a hyperkahler manifold acted on by SU(2) (the 
action descends to 5*0(3) for AH). The manifolds have a 'core region' and an 'asymptotic 
region' with an asymptotic action of 0(2). In addition to the 2-parameter family of complex 
structures given by the hyperkahler metric (and rotated by SU(2)) there is another complex 
structure compatible with the SU(2) action, but giving the opposite orientation to that of 
the 2-parameter family. For TN this complex manifold is just C 2 and for AH it is CP 2 \IRP 2 . 

In both cases the asymptotic 0(2) action gives us topologically a disc bundle at infinity and 
this can be identified with the normal bundle of the surface X that naturally compactifies 
the manifold. For TN this compactification is CP 2 with X = CP 1 , and for AH it is CP 2 with 
X = MP 2 . We have shown that these give rise to opposite signs: 

For TN X 2 = +1 while for AH X 2 = -1 . 

As explained in the main text, TN and AH are our geometric models for the electron and 
the proton, and we identify the self-intersection number of X with minus the electric charge. 
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B Metric properties of TN and AH 



B.l Coordinates and conventions 

In this paper we need to parametrise SU (2) explicitly at various points, and we use generators 

* i= ^ ti= K< o)' t2= ^ r2= K-i o)' t3= ^ r3= Ko -<)' (ri) 

where r , a = 1,2,3, are the Pauli matrices; the commutators are [t a ,tf,] = —e a b c t c . We use 
Euler angles 9 G [0,7r), G [0, 27r) and -0 G [0, 47r) for parametrising 577(2) as follows: 

i_ e2 vv w sin ^ e 2 cos lq i 
Defining left-invariant 1-forms r\ a on SU (2) via 

g - l dg = rati + r} 2 t 2 + r} 3 t 3 , (B.3) 

we compute to find 

rji = — sin ip dO + cos ip sin # d<p , 
?72 = cos ip d9 + sin ^ sin # <i0 , 

?7 3 = dip + cos 6> # , (B.4) 

satisfying drji = \(-ijkf]j A r] k . They descend to 5*0(3) by simply restricting ip to [0, 2n). 

Explicitly, generators of the Lie algebra of 5*0(3) are (T a ) bc = t ahc . They also satisfy 
[T a , 7&] = —e a b c T c . We can then parametrise 50(3) matrices via 

i?(0,^) = e^e e7 V T3 , (B.5) 

with 6 G [0, tt), G [0, 2tt) and ^ G [0, 2tt). Then 

R^dR = + V2 T 2 + m T 3 (B.6) 

provides an alternative definition of the left-invariant 1-forms on 50(3): one obtains the 
same expressions as in (B.4), but with the range of angles automatically appropriate to 
50(3). 

Both TN and AH can be parametrised in terms of Euler angles and a radial coordinate. 
For TN, the angular ranges are G [0, 7r), <fi G [0, 2n) and ip G [0,47r). For AH they are 
6 G [0, 7r), (p G [0, 2n) and tp G [0, 2n) with the additional Z 2 identification 

(8,<i>,rl>)~(n-8,<i> + n,-il>), (B.7) 
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which, in the asymptotic region, is the simultaneous reversal of spatial and fibre direction. 
The following angular integrals, which enter into various calculations in this paper, are 
therefore 



r] 1 Ar] 2 Ar] 3 = -16n 2 , 



SU(2) 

f 

50(3) 



f]\ A 7] 2 A 7] 3 = -87T 2 , 



/ 77! A 772 A % = -4vr 2 . (B.8) 

J(SO(3)/Z 2 ) 

Note that 771 is invariant under (B.7), but 772 and 773 change sign. 
Recalling that the metrics on TN and AH are of the Bianchi IX form 

ds 2 = firfdr 2 + airfril + b(r) 2 r]l + c{r) 2 r, 2 3 , (B.9) 

we introduce the tetrad 

e 1 = a77x , 9 2 = br] 2 , 9 3 = C773 , 9 A = fdr . (B.10) 

The self-duality of the Riemann tensor computed from (B.9) with respect to the volume 
element 

dV = ar/i A br] 2 A cr] 3 A fdr = —fabcdr A 771 A 772 A 773 (B-H) 
is then equivalent to the set of ordinary differential equations 

= (b - c) 2 - a 2 + 2\bc , + cyclic, (B.12) 

j dr 

where '+ cyclic' means we add the two further equations obtained by cyclic permutation of 
a, 6, c, and A is a parameter which has to be either or 1. In all cases the resulting metrics 
are hyperkahler, but in order to obtain metrics whose hyperkahler structures are rotated by 
the SU (2) action we need to set A = 0. This is the case for both TN and AH, so A = in 
(B.12). 

One checks the equivalence of the self-duality of the Riemann tensor and (B.12) as follows. 
Denoting the Riemann tensor by R, its component 2-forms with respect to (B.10) are 

R12 = d/J> 3 A 773 + (/i 3 - /ii/^2 - AiA 2 )t7i A 772 , 

R u = d\ 3 A 773 + (A 3 - Ai/i 2 - A 2 /!i)t7i A 772 , (B.13) 
and similar expressions for the components obtained by cyclic permutation of 1, 2, 3, where 

b 2 + c 2 -a 2 c 2 + a 2 -b 2 a 2 + b 2 -c 2 

A*l = 7jZ > ^2= o ' ^ 3 = o~k B - 14 

2oc 2ca 2ab 
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and 

^ 1 da ^ _ 1 db ^ 1 dc 
1 / dr ' 2 / dr ' 3 / dr 

When (B.12) hold, then Aj = //j + A — 1, i — 1, 2, 3, with A = or A = 1, and the self-duality 
relations 

R12 = -R34 , -R23 = Ru , R31 — R24 (B.16) 

are easily verified. One also checks that each of the three independent 2- forms R12, R23, R31 
is self-dual with respect to (B.ll), as expected. 

The only solutions of (B.12) with A = which give rise to complete manifolds whose generic 
577(2)- or 50(3)-orbit is three-dimensional are the TN and AH metrics, whose coefficient 
functions we discussed in Sect. 3. It is important that the coefficient function c in (3.1) is 
negative for all r in the AH manifold. It implies in particular that the canonical volume 
element/orientation dV in (B.ll), with / = — b/r, 

ab^c 

dV = -fabc dr A m A m A m = sin 9 dr A dO A # A dib , (B.17) 

r 

has opposite signs for TN and AH, in the following sense. Assuming the same orientation 
Vi A^2 A773 of the 577(2) and 50(3) orbits in TN and, respectively, AH, the radial direction is 
oppositely oriented in the two cases: the natural radial line element —fabcdr has the same 
orientation as dr for TN, but the opposite orientation for AH. This will be important when 
evaluating various integrals. It means that, when using the coordinates r,9,(p,ip we can 
use the conventional orientations for these coordinates when computing for AH, but should 
use the opposite orientation when computing for TN. Thus we integrate in the negative 
r-direction when calculating on TN. 



B.2 L 2 -norms of the Riemann curvature 

The squared L 2 -norm of the Riemann tensor is 

\\R\\ 2 M = [ 2(it42 A R12 + R23 A -R23 + R31 A -R31) , (B.18) 

where M is TN or AH, and we have used the self-duality of the curvature. In terms of 
the functions (B.14) the integrand can be simplified (see also [38], where this calculation is 
carried out for the AH metric) and becomes 

2{R l2 A R 12 + R23 A R 23 + R31 A R31) = dFAr ]l Ar ]2 A 773, (B.19) 

with 

F = 2(/ix + 1^2 + 1^3 — I) 2 — 8 / u 1/ u 2 /i2 • (B.20) 
We denote these functions in the TN and AH cases by F TN and Fah- 
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For TN, we can compute explicitly and find 

^TN(r) = -8^±| ¥ . (B.21) 

Integrating (B.19) and using (B.8) as well as our sign convention for the volume form we 
find 

||i?|| 2 N = (-16tt 2 )(F tn (0) -Ftn(oo)) = 8tt 2 . (B.22) 

For AH, similarly, 

II^IIah = (-47r 2 )(F AH (oo) - F ah (tt)) . (B.23) 

To complete the calculation we need the behaviour of Fah at infinity and at ir. One finds 
for AH 

2 

//i(oo) =// 2 (oo) = ^(oo) =0, ^3(00) = 1 - ^2(00) = 1, (B.24) 

as well as 

1 1 

Hi(ir) = -1, / i 2 ( 7 r) = -, /i 3 (7r) = -, (B.25) 
where we took careful limits, using (3.7) and (3.8). Therefore 

F AH (oo)=0, F AH (vr)=4, (B.26) 

so 

I \R\ \l R = (-4tt 2 ) x (0 - 4) = 16tt 2 . (B.27) 

C Harmonic forms on TN and AH 
C.l Rotationally symmetric harmonic forms 

The question of computing harmonic 2-forms on the TN and AH manifolds arose in physics 
in the context of testing S-duality, see [33] for AH and [30, 31] for TN. Harmonic forms also 
play a role as curvatures of line bundles over TN and AH, see [32] for a discussion of this for 
AH where the form later used for testing S-duality [33] is the curvature of an index bundle. 
This appendix contains a systematic discussion of these forms and also a harmonic form on 
the branched cover of AH which has not previously been considered in the literature. The 
physical interpretation of the forms as electric and baryonic fluxes is given in Sect. 8. 

On both TN and AH, rotationally symmetric 2-forms can be written in terms of the metric 
coefficient functions /, a, b, c appearing in the respective metrics and functions f\, fi and f 3 
of the radial coordinate r as 

uf = ff (br] 2 A o] 3 ± arji A fdr) , 
w 2 = ft ( cr ?3 A arn ± br/2 A fdr) , 

= jf (arj! A br/2 ± cr/ 3 A fdr) , (C.l) 



33 



with + standing for self-dual and — standing for anti-self-dual. Introducing the functions 

gf = f?bc, g± = f±ca, gf = f±ab, (C.2) 
closure of the forms (C.l) implies the ordinary differential equations 

dgf af ± 



dr 6c 

^2 &/ 



± 



T <?2 , 

ar ca 

For solutions of these differential equations (we shall discuss below in which cases regular 
solutions exist), the forms (C.l) can be written 

vf = gfv2 A r/ 3 + dgrf A 77! , 
^ = 02^3 A ?7i + dffa A rj 2 , 

^ = #3^1 A V2 + ^3 A 773 • (C.4) 
All these forms are locally exact, that is, 

wf = disfH) , * = 1,2,3, (C.5) 

but the 1-forms in brackets may not be globally defined and, when defined, they may not be 
L 2 , so that the corresponding forms uf are not necessarily L 2 -exact. 

We are interested in regular and bounded solutions of the equations (C.3). Using our 
convention / = —b/r, the functions appearing in the differential equations are then 

af a bf _ b 2 cf _ c 

— — , — , — — . l^-DJ 

be rc ca rca ab ra 

C.2 Taub-NUT 

Substituting the TN coefficient functions (3.5) (with e = 1 and m — 2), the functions (C.6) 
simplify further to 

_^_ = _*!_ = _^ c - 2 (C.7) 
rc rca 2r ' ra r(r + 2) 

It is easy to check, and was shown in [30] and [31], and earlier in [28] and [29], that only the 
equation for g£ has a solution which is both regular at the origin and bounded as r — > 00. 
Integrating, one finds explicitly 

9t(r) = C^, (C.8) 
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with an arbitrary constant C. Setting C — 1, (C.4) gives the associated normalisable, 
harmonic form 



Then, using (B.8), 



u t A u>t = / — — rfr x (-167T 2 ) = 16tt 2 . (C.IO) 

TN </oo V + A ) 



C.3 Atiyah-Hitchin 



On AH, the analysis of harmonic forms is analogous. We write Gf for the solutions of (C.3) 
with the radial functions a, b, c, / of the AH metric, and write the forms obtained by solving 
these equations as 

Qf = Gfr} 2 Am + dGf A 771 , 
nf = G^Vs A r]i + A T] 2 , 

= Gfr ]l Ar ]2 + dGf A V3 . (C.ll) 

As for TN, all these forms are formally exact, 

nt = d(Gt Vl ), i = 1,2,3, (C.12) 

but the 1-forms in brackets may not be globally defined, as we shall see. 

Note that only the forms Qf respect the identification (B.7). However, as explained in the 
main text we also consider the branched double cover AH of AH and therefore keep all forms 
in the discussion. On AH, we use (3.8) and (3.7) to find that the functions appearing in the 
differential equations (C.3) have the following behaviour near the core (i.e. (r — ir) small), 

rc 7i z rca ra 2{r — ir) 

and the following behaviour for large r, 

a b 2 1 / 2\ c 



rc rca 



1 1 — — J , , -. (C.14) 

2 V r J ra r(r - 2) y J 



For the differential equations (C.3) on AH, this means that only the equations for G^ and 
have solutions which are finite at the core and remain bounded for large r. Both of these 
solutions are discussed and interpreted in the main text. 



For any solution of (C.3) on AH, we note 

/ n ± Ail ± = [ rf((G ± ) 2 ) A ^ A 772 A 773 
J ah Jar 



;-4tt 2 ) ((G±) 2 {oo) - (G±) 2 (*)) , (C.15) 
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where and G ± stand for any of the forms and coefficient functions in (C.ll), and we 
again used (B.8). We also note that the L 2 -norm of each of these forms can be written as 



|f2 ± || 2 = ± / ft* A ft* 



pdV, (C.16) 



AH 

where dV is the volume element defined in (B.ll) and p is a density which we interpret as 
an energy density in the main text. Its general form is 

which can be simplified in each of the cases, using the differential equations (C.3) on AH. 
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